ENERGY LOSS OF HIGH p ± HADRONS 
BY FINAL HADRONIC STATE 



in 
o 



K. GALLMEISTER, C. GREINER AND Z. XU 

Institut fiir Theoretische Physik, Universitat Giessen 
~^\> Heinrich- Buff- Ring 16, D-35392 Giessen, Germany 

o 
o 

CN 1 Abstract 

F ^ \ In the context of the 'jet quenching' phenomena typically materializa- 

pLn ' tion of the jet is assumed to take place in vacuum outside the reac- 

■ tion zone. On the other hand quantum mechanical estimates give a 

hadronization time on the order of only a few fm/c for jets materializ- 
ing into hadrons with transverse momenta of p± < WGeV, which thus 
should well take place inside the fireball. Typical (in-)elastic collisions 
of these high pj_ particles with the bulk of hadrons of the fireball have 
a rather low invariant mass and are thus nonperturbative. An analysis 
£N| ' within an opacity expansion in the number of collisions by means of 

the FRITIOF collision scheme for various hadrons will be presented. 
\ It shows that late hadronic collisions can substantially account for the 

modification of the high pj_-spectrum observed for central collisions at 
RHIC. 

1 Introduction: Qualitative estimates 

As one of the very interesting first results, RHIC experiments have established 
a significant suppression of moderately high p± hadrons produced in central 
A+A collisions compared to rescaled peripheral collisions or rescaled (and ex- 
trapolated) p+p collisions [|TJ . The results are a clear hint for nuclear medium 
effect (s) at work. 

The most popular explanation is the onset of the occurrence of so called 
'jet quenching' f2|. It is the idea that a high energetic parton moving through 
a color dense (and then deconfined) medium will lose considerable energy most 
dominantly by gluon radiation, and so its final fragmentation will give rise to 
a leading hadron with considerable lower energy. A recent and much involved 
calculation based on the GLV formalism for a finite number of colored collisions 
(and including the Cronin effect as well as a slight modification of the gluon 
distribution due to shadowing) provide a good agreement with the data. Still, 
though, a phenomenological opacity parameter has to be adjusted [El. It has 



also been suggested that simple incoherent scattering in the further partonic 
evolution of the fireball should be of relevance |4| . Contrary, also the proposal 
was made that the observed spectra for central collisions show a significant 
m± scaling being a manifestation of a direct remnant of the initial gluons 
which were liberated from gluon saturated nuclear distribution functions ||. 
Here, any possible final state interactions which might or should alter the 
distributions are completly discarded. 

A major concern, and which lead to the present study, is that it is typically 
assumed in all these descriptions that the partons somehow exits the collision 
region before finally fragmenting into a 'jet' of hadrons. The potential mag- 
nitude of the hadronization time (or, to be more precise, the time to build 
up the hadronic wavefunction) is based on a relativistic and simple quantum 
mechanical estimate for (a) light or (b) heavy quarks |J 

(a) t^ adr « ER 2 <-> (b) t h a adr « R . 
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Taking (a) with the average radius of the pion R n « 0.5 fm or (b) substituting 
tuq by nip for a p-meson and taking R p rs 0.8 fm, one has for the formation 
time a crude understanding given as 

t F « 1 - 1.2 (E/GeV) * fm/c . (1) 

Hence, for leading hadrons with moderately high p± < lOGeV original point- 
like jet-partons have established already a complete nonperturbative, transver- 
sal wavefunction after traveling a distance in length of smaller or equal than 
10 fm. Accordingly, the jets should, to a large fraction, materialize into 
hadrons still inside the expanding fireball, which has a transversal size of 
R ~ 8/m + 0.5 * ct, where t denotes the local proper time of the expand- 
ing system. 

The state becoming the leading hadron will interact as a (pre-)meson with 
the surrounding low momentum mesons with a cross section of roughly a(t) = 
(To * (t/tp) as long as it is not completely being build up ||. cr ~ 10 — 15m6 
denotes the total cross section of two mesons. The density of hadrons in the 
late fireball changes from about 1 fm~ 3 to 0.1 fm -3 . The mean free path of the 
fast hadron is then estimated to be A ~ 1 — 10 fm, meaning a few collisions 
L/X — 0, 1, 2, 3 — The system is potentially rather opaque! 

Looking at the available energy for an individual collision (compare Fig. 1), 
even for a value of p± = 10 GeV, one gets a y/s < 2 GeV, if the target is a 
pion at rest, and only with a p as target, one has an invariant mass of above 
2.5 GeV for particles with a transverse momentum larger than 3 GeV. 
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Fig. 1: The CMS energy for three targets at rest (ir,K,p, depicted by the labels near the 
curves) and three types of projectiles (n, K, n/p, indicated by the solid, dashed and long 
dashed lines) as function of the momentum of the projectile. 

Hence, for the considered transverse momentum region, one has either elas- 
tic scattering, resonance scattering or also inelastic scattering resulting in a few 
final hadrons. In the following we treat the collision with the FRITIOF Monte- 
Carlo scheme with a p-meson as characteristic target hadron being at rest. To 
also stress some model dependence, we later compare with elastic scattering 
on either a p or a ir as target particles, where the scattering is taken simply as 
isotropic f7|. 

2 Energy loss by (multiple) final state hadronic inter- 
actions 

The energy loss for one single interaction of a leading hadron i with momentum 
p± on a low transverse momentum target (taken as ~ 0, and both hadrons 
being located at the same space-time rapidity and momentum rapidity in a 
Bjorken-type picture), going into a hadron j with (lower) momentum p± can 
be written as 

fjiPl-) = 2 / d P±° fi(Pu) 9ij(PUiP±) , ( 2 ) 

i 

while fi(p± ) is the initial distribution, gij(p± ,p±) is a folding matrix and 
fj(p±) is the resulting distribution. i,j stand for the particles n°, n ± , K°, K° , 
p, p, n and n. The folding matrix depends on the model employed for 
the collision (for further details we refer to 0). 
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Fig. 2: Comparison of the calculations at midrapidity (for one unit in rapidity) with the 
(preliminary) PHENIX data |j| on charged hadrons for three different centralities. The thin 
dashed lines indicate the results from PYTHIA, scaled with a corresponding binary collision 
number N co u. The thick solid line depicts the rate for the most central region where the 
particles have suffered on average one hadronic collision. 

For the initial distribution ff of individual hadrons i at moderately high 
p± when entering the final hadronic stage we have no direct information. All 
various ideas, as mentioned at the beginning, might contribute and already 
steepen the spectrum. In order to see the possible effect of the late hadronic 
interactions, we use individual distributions generated by PYTHIA v6.2 |5J. 
The parameters were first adjusted to UAl-results of pp-collisions. Isospin 
averaged N-N-collision were then calculated (for details we refer to 0). With 
the appropriate scaling of binary collisions (N co u = 2.4 for 80. . . 92%, 12.6 
for 60. . . 80% and 945 for 0. . . 5% centrality), the calculations are in perfect 
agreement with the preliminary PHENIX data for charged hadrons and for 
neutral pions in the case of peripheral collisions, as can be seen from Figs 2 
and 3. 

Moreover, one can see a dramatic effect when employing on average one 
single final inelastic collision modelled via the FRITIOF scheme. In case of the 
charged hadrons, one single collision already would slightly overestimate the 
suppression. On the other hand one collision is just right to perfectly explain 
quantitatively the modification of the momentum distribution of the 7r ! (We 
have only shown the modifications for p± > 2 GeV, as the spectrum below 
should be dominated by soft physics and hydrodynamical expansion.) 
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Fig. 3: Like Fig. 2, but now a comparison of the present calculations with the PHENIX data 

U on 7T . 

Hence, any potential energy loss by final state hadronic reactions is in 
the same range as considered for the jets in deconfined matter. In Fig. 4 we 
show calculations for various number of collisions now at a ^/snn = 200 GeV 
up to p± = 10 GeV. The three employed models (for treating the collisions) 
do give more or less similar predictions. The spectrum steepens considerably 
for a larger effective opacity < L/X >. On the other hand at some stage 
the complete evolution of all hadronic particles has to be incorporated to see 
in full detail the dynamical interplay among 'soft' and 'hard' hadrons. The 
spectrum can not become more steepen than exponential as then the system 
has thermalized also in the high p± hadronic degrees of freedom. 

In summary, we have motivated that (pre-)hadrons stemming from a jet 
should still materialize in the dense system for momenta up to 10 GeV. The 
late hadronic interactions have a clear and nonvanishing effect in suppressing 
the spectrum. On the average one single such interaction should already be 
enough to explain quantitatively the RHIC results. Our finding, however, sig- 
nals a warning which has to be addressed in future studies very accurately. 
'Omnes viae Romam ducunt': Various different roads at present can account 
for the possible explanation of the modification in the spectra for central col- 
lisions. Deductions for possible QCD effects of a deconfined QGP phase on 
the materializing jets have to disentangle from the here investigated final state 
interactions, before definite conclusions on the importance of a potential dense 
partonic phase, or any other effect, can quantitatively be drawn. 
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Fig. 4: Resulting p± -spectrum of charged hadrons at midrapidity for ^snn = 200 GeV and 
for < L/X >= 0, 1, 2, 3 collisions. Hadronic collisions are treated in three different ways 
(compare text). 
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